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We report on and characterize, via molecular dynamics (MD) studies, the evolution of the struc- 
ture of CusoZrso and Cu64Zr36 metallic glasses (MGs) as temperature is varied. Interestingly, a 
percolating icosahedral network appears in the Cue4Zr36 system as it is supercooled. This leads 
us to introduce a static length scale, which grows dramatically as this three dimensional system 
approaches the glass transition. Amidst interpenetrating connections, non-interpenetrating connec- 
tions between icosahedra are shown to become prevalent upon supercooling and to greatly enhance 
the connectivity of the MG's icosahedral network. Additionally, we characterize the chemical com- 
positions of the icosahedral networks and their components. These findings demonstrate the im- 
portance of non-interpenetrating connections for facilitating extensive structural networks in Cu-Zr 
MGs, which in turn drive dynamical slowing in these materials. 



I. INTRODUCTION 

Understanding the structure-property relationships of 
metallic glasses (MG) is a pressing matter in the field 
of material science 1 - - — . The culmination of past research 
reveals that the vast diversity of MG species, and their 
properties, are rooted in an equally-broad set of struc- 
tural archetypes 6 -"— . Perhaps even more cumbersome is a 
decades old outstanding problem in the study of glasses- 
namely explaining the extreme slowing of dynamics in a 
liquid as it is supercooled towards a glass transition^—. 
Many theories, such as the theory of random first order 
transition a 10 ' 11 , spin glass approaches^, and others^— 
predict the existence of rapidly growing length scales 
(in particular, a static length scale) to accompany the 
marked increase of the viscosity of supercooled liquids. In 
fact, there are rigorous results predicting the appearance 
of such lengths^ 3 -. While much progress^— has been 
made in studying the structures of various glass-forming 
systems, to date, no notable increase in standard static 
length scales has been observed in most studies. 

Cu-Zr is a popular representative of Early Transition 
Metal-Late Transition Metal MGs; this is, in part, be- 
cause it has a high glass forming ability (GFA) for a 
broad range of compositions^—. The binary compo- 
sition of Cu-Zr reduces the complexity of the possible 
local atomic structures, making this system ideal for the 
study of the evolution of the spatial structure of liquids 
as they are supercooled to form a glass. Recent research 
efforts demonstrate a structural hierarchy within Cu-Zr 
MGs that appears to be central to its structure-property 
relationships^^—. Strong evidence has been found to 
suggest that the presence of Cu-centered full icosahedra 
is uniquely responsible for dynamical slowing during the 
formation of the Cu-Zr MGs^^^. Aside from their full 
and distorted icosahedral Voronoi signatures, the Voronoi 
polyhedron landscapes of Cu-Zr in the liquid and glass 
phases resemble one another considerably 4 ^. The icosa- 
hedral clusters often interpenetrate one another, so that 
five atoms coincide on two icosahedral-shells. These in- 
terpenetrating connections of icosahedra (ICOI) exhibit 



strong spatial correlation a 34 ' 47 and produce stable string- 
like networks of icosahedr a 34 ' 46 " —. Extensive ICOI clus- 
ters are found to posses a high average elastic rigidity^ 6 - 
and to enhance more general mechanical properties of 
Cu-Zri2 MG. 

In light of the importance of icosahedral clusters and 
ICOI in shaping Cu-Zr MG, it is pressing to consider 
the roles of non-interpenetrating connections in order to 
obtain a complete picture of icosahedral networking in 
this amorphous alloy. To date, the prominence of these 
connections has only been noted 3 ^. In this paper, MD 
simulations are used to study non-interpenetrating con- 
nections, and to characterize their effect on the connec- 
tivity and composition of icosahedral networks in Cu-Zr 
MGs, using CusrjZrso and Cu64Zr36 as models. Our anal- 
ysis reveals that, near the glass transition, a large num- 
ber of non-interpenetrating connections develop amongst 
ICOI-structures and significantly enhance the unification 
of icosahedral networks in these MGs, thus impacting 
their mechanical properties. Considering the full icosa- 
hedral network affords us the rather unique opportunity 
to identify, in Cu64Zr36, a percolating icosahedral network 
that threads the entire glassy system, providing a large 
scale static structure which grows rapidly as the system 
approaches the glass transition. We also suggest that the 
non-interpenetrating connections are particularly influ- 
ential in MGs like CusoZrso, which exhibits relatively few 
ICOI and yet is a good glass former. Moreover, with all 
of the connection types being included, we reveal that the 
icosahedral network and the otherwise liquidlike matrix^ 
have significantly different chemical compositions. For 
example, vertex and face-sharing connections preferen- 
tially incorporate Zr atoms, increasing the density of the 
icosahedral network. For both alloys, the icosahedral net- 
work is Zr-rich, relative to the overall sytem composition, 
and is immersed in a Cu-rich liquidlike matrix. 

In section II, we will proceed to discuss the methods 
for our models of Cu-Zr, MD simulations, and structural 
analysis. Section III characterizes the basic icosahedral 
ordering in our systems as they are cooled from the liq- 
uid phase to the glass phase. In section IV we mea- 
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sure the population and distribution of connections in 
the icosahedral network. Sections V, VI, and VII, re- 
spectively, presents weighted connectivities, introduces a 
static length scale, and details the chemical compositions 
of the icosahedral networks for CusoZrso and Cu64Zi - 36. 



II. MOLECULAR DYNAMICS SIMULATIONS 
AND VORONOI TESSELLATION 

The MD simulations conducted in this study use a 
semi-empirical potential developed by Mcndclcv ct al, 
which describes many-body interactions using the Finnis 
and Sinclair (FS) generalization of the embedded atom 
method (EAM)£i. It was created, in part, by being fit- 
ted to x-ray diffraction data, liquid density, enthalpy of 
mixing data, and first-principles results. The potential 
was developed with an emphasis on describing both liq- 
uid and amorphous Cu-Zr alloys, and was shown to pro- 
vide a very good description of the structures and some 
properties of the liquid and glass phases of Cu64.sZr35.s. 
An earlier version of this potential^ was also shown to 
produce reliable atomic structures for Cu64.5Zr 35 . 5 MG 
when compared to atomic structures generated by the 
reverse Monte Carlo (RMC) method, which was con- 
strained by independent x-ray diffraction data and ab 
initio MD simulation results^. The semi-empirical po- 
tential used in this paper accurately recreates the prop- 
erties of pure Zirconium and is expected to yield reason- 
able structures for CusoZrso as well. Accordingly, the 
structure factor data and Voronoi polyhedra distribution 
obtained for this composition are in good agreement with 
the published x-ray and Voronoi polyhedra data of M. Li 
et a/4£(not shown). 

Canonical (NVT) systems of CusoZrso and Cu64Zr36 
were simulated using LAMMPS^. Each simulation con- 
sists of 2000 atoms contained in a box with periodic 
boundaries. The initial configuration of an alloy is gen- 
erated randomly, and the alloy is allowed to melt and 
evolve for .25 ns at 1600K, using 5 fs time steps. The al- 
loy is subsequently quenched to its target temperature at 
a rate of 10 12 K/s, with the volume of the simulation box 
set to yield zero pressure for this temperature. Finally, 
the quenched alloy evolves for 1 ns at the target temper- 
ature, during which 78 snapshots of the atomic configu- 
ration are recorded. The target temperatures used range 
from 300 K to 1600 K in steps of 100 K. The process is 
repeated so that, ultimately, three simulations with in- 
dependent initial configurations were carried out for each 
of the fourteen target temperatures. To check for system 
size effects, canonical systems with 10,000 atoms and a 
target temperature of 300 K were also simulated for both 
alloy compositions. The approximate values of the glass 
transition temperatures (T g ) for CusoZrso and Cu@4Zr36 
were found to be 700 K and 800 K, respectively. 

An analogous NPT simulation (P = 0) of Cu 50 Zr5o 
with a target temperature of 400 K was conducted for 
the sake of comparison. The resulting average partial 
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FIG. 1. (Color Online) Average percentage of atoms in a 
Cu-Zr alloy that participate in a 13-atom icosahedral cluster 
(center atom included). For T = 300 K, data points are in- 
cluded for a system size of N — 10, 000 atoms in addition 
to N — 2, 000 atoms. The error bars indicate the standard 
deviation of the value. The glass transition temperatures are 
marked for the respective compositions. 

pair correlation functions (ppes) were found to be sta- 
tistically identical to those from the corresponding NVT 
simulations. 

The atomic structures in the MD-generated alloys were 
characterized using Voronoi polyhedron-analysis for each 
atoro^ i 54 ' 55 . Using a procedure similar to that of Hsu 
and Rahman, cluster vertices that were separated by dis- 
tances less than seven orders of magnitude of the simula- 
tion box edge length were consolidated into one vertex^. 



III. ICOSAHEDRAL ORDERING IN Cu-Zr 
ALLOYS 

To study the development of icosahedral ordering in 
Cu-Zr, we consider the results of our Voronoi tessella- 
tion analysis: a Cu-atom with a Voronoi index (0,0,12,0) 
and its 12 Voronoi neighbors form a full icosahedral clus- 
ter. The index (71.3, rc.4, n&, ...) specifies the number of 
faces, rife, with k edges present on the Voronoi polyhe- 
dron. Zr-centered full icosahedra are seldom observed 
amongst the Voronoi polyhedra in these alloys, which is 
consistent with previous MD studies, thus (0,0,12,0) clus- 
ters will henceforth refer only to Cu-centercd full icosa- 
hedra. Figure [T] displays for a range of temperatures the 
percent-abundance of atoms involved in (0,0,12,0) clus- 
ters. In their liquid phases above 1300 K, fewer than 10% 
of the atoms participate in (0,0,12,0) clusters for both 
compositions, but supercooling yields a great enhance- 
ment in the degree of icosahedral ordering. Accordingly, 
roughly 28% and 53% of the atoms of glassy Cu 50 Zr 5 o 
and Cu64Zr3g, respectively, participate in the structure 
of full icosahedra (including their center Cu-atoms) . This 
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FIG. 2. (Color Online) [(a) and (b)] The average percentage of (0,0,12,0) clusters involved in: at least one ICOI connection, 
strictly non-ICOI connections, no connections for (a) Cu64Zr36 and (b) CusoZrso. The black symbols correspond to a system 
size of N = 10, 000 atoms, [(c) and (d)] The average percent distribution of connection types between icosahedra in (c) Cu^Ztsq 
and (d) CusoZrso. The cluster snapshots in the right column of the figure depict the different types of connections between 
icosahedra. The red spheres represent the center Cu-atoms of the two icosahedra involed. 



compositional trend of an increase in icosahcdxal order- 
ing, and subsequently an increase in T g , that accompanies 
an increase in the fractional abundance of Cu-atoms has 
been demonstrated and studied previousl y 40 i 41 . 

IV. ATOMIC CONNECTIONS BETWEEN 
ICOSAHEDRA 

Below T g , we observe that approximately 8.6% of 
Cii64Zr3e's Voronoi polyhedra in our simulation snap- 
shots are (0,0,12,0) clusters. If these clusters were iso- 
lated from one another, it would require that 112% of the 
system's atoms participate in their structure. Consider- 
ing that only about 53% of this alloy's atoms are actually 
involved in these clusters, it is obvious that many clusters 
share atoms amongst one another and are thus joined to- 
gether. Therefore, two (0,0,12,0) clusters are defined to 
be 'connected' if they share at least one atom. The inclu- 
siveness of this definition will be justified as we further 
study the details of icosahcdral structures. Icosahedra in 
our simulated Cu-Zr alloys are observed sharing vertices, 
edges, and triangular faces in addition to exhibiting the 
interpenetrating connections described earlier. Examples 
of these connections are depicted in Figure [H We will 
study the creation of these connections between icosahe- 
dra during the process of quenching to understand how 
they form the mechanically-stiff icosahedral networks. 

Figure [2] quantifies the development of connections be- 
tween icosahedra as the Cu-Zr alloys arc quenched from 
the liquid phase. Figures^ and [2b illustrate that as the 
number of icosahedra in a Cu-Zr alloy increases rapidly 



during supercooling (refer to Figure 1), the proportion 
of those participating in these stable ICOI grow compa- 
rably. Above 1300K, less than a third of the icosahedra 
are involved in ICOI for both alloys. The prevalence 
of ICOI-based structures is heightened dramatically dur- 
ing supercooling, and once Cu64Zr36 and CusoZrso are 
arrested in their glass phases, approximately 80% and 
60%, respectively, of icosahedra have at least one inter- 
penetrating connection. Very few icosahedra are isolated 
in glassy Cu-Zr MGs, e.g, about 1.7% and 14%, respec- 
tively, in Cu 64 Zr 36 (Fig. 2a) and Cu 5 oZr 5 o (Fig-UJa). Fig- 
ures and [2Ji display the relative abundance of each 
connection- type in the alloys. They show that vertex 
and face-sharing connections become increasingly preva- 
lent as the alloys are quenched and enter their glass 
phases. For example, below T g in Cug4Zr36, these con- 
nections are as common are ICOI (Fig. &p). This is in 
stark contrast to the aforementioned results of Figures^ 
and [2]3, indicating that, as the density of icosahcdral 
clusters and the number of ICOI increase during super- 
cooling, vertex and face-sharing become prevalent amidst 
ICOI-based structures. We will further study how these 
non-interpenctrating connections are incorporated into 
the icosahedral structures of Cu-Zr MG alloys. 

Above 1300 K, the connections between icosahedra are 
constantly forming and breaking (as are the icosahedra 
themselves) so that no single connection between two 
icosahedra is observed in consecutive snapshots (sepa- 
rated by 13 ps) of a simulation. Viewing snapshots of 
the liquid alloys and assessing the average abundances of 
the different connection-types thus provides a measure of 
the relative prominences and stabilities of their resulting 
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FIG. 3. (Color Online) [(a) and (b)] The average distribution of the number of interpenetrating connections per icosahedron 
in (a) Cu64Zr36 and (b) CusoZrso MGs at 300 K (N — 10, 000 atoms). The error bars indicate the estimated error in the mean 
value, [(c) and (d)] The average distribution of the number of connections (vertex, edge, face-sharing, or interpenetrating) per 
icosahedron in (c) Cue4Zr36 and (d) CusoZrso MGs at 300 K (N = 10,000 atoms). These are stacked-histogram plots in which 
the blue bars correspond to icosahedra that participate in at least one interpenetrating connection, and the red bars correspond 
to those that do not participate in ICOI. 



structures. Figures [2}; and [2]i show, for temperatures 
above 1000 K, that ICOI comprise more than 40%, a 
significant majority, of the connections present in both 
alloys. This provides strong evidence that the ICOI is 
in fact more robust than the other types of connections 
and leads us to consider its role in the liquid alloys as 
they cool. It is incorrect, however, to think of the glassy 
icosahedra! network as simple chains of interpenetrating 
icosahedra, as our considerations of Figures [5J: and [2ji 
below T g demonstrate that these ICOI structures also 
incorporate a significant number of non-interpenetrating 
connections, as we elaborate on next. 



V. NON-INTERPENETRATING 
CONNECTIONS AND NETWORK 
CONNECTIVITY 

To study the structures of the icosahedral networks 
found in Cu-Zr MGs, we begin by characterizing the for- 
mative ICOI-structures. Figures [3h and [5}} show the av- 
erage distribution of the number interpenetrating connec- 
tions per icosahedron in the alloys at 300 K. As expected 
from our considerations of Figure [2j Cu64Zr36 contains 
more extensively connected ICOI-structures than does 
Cu5oZi'5o MG with an average of 2.2 interpenetrations 
per connectcd-icosahedron versus 1.5, respectively. Both 



compositions primarily exhibit chain-like and triangular 
arrangements of interpenetrating icosahedra. These con- 
nectivity measurements are corroborated and elaborated 
on by earlier studie a 42 ' 47 . Figures and [3Ji incorpo- 
rate all connection types, and provide a uniquely detailed 
characterization of the icosahedral-structures found in 
these MGs. They display the average distribution of 
the total number of connections per icosahedron, while 
distinguishing icosahedra that participate in ICOI from 
those that do not. From this perspective, the average 
numbers of connections per interpenetrating icosahedron 
increase to 6.8 and 3.0, respectively, for Cu64Zr36 and 
CusoZrso. These distributions thus describe more intri- 
cate, highly connected networks than do to those in Fig- 
ures [3^ and Furthermore, it becomes clear that in- 
terpenetrating icosahedra frequently participate in other 
connection types and that they serve as the nodes with 
the highest connectivity in the icosahedral networks of 
the MGs. For example, in Cu64Zr36 a sizable fraction 
of ICOI-icosahedra can be found participating in ten 
or more connections, while non-ICOI-icosahedra are sel- 
dom found with more that six connections. It follows 
that significant number of non-interpenetrating connec- 
tions must occur within and between ICOI-structures. 
These highly connected nodes in the networks likely cor- 
respond to icosahedra with very low atomic mobilities, 
which arc known to be responsible for dynamical slow- 
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FIG. 4. (Color Online) Simulation snapshots of icosahedral 
structures containing a) pinning and b) bridging connections. 
Only the central Cu atoms of the icosahedra are displayed. 
Black lines are drawn between center Cu atoms to indicate 
non-ICOI (vertex, edge, or face-sharing) connections between 
their respective icosahedra. 

ing in these MG o 34 ' 41 ' 42 . When inspecting such densely 
packed icosahedra, simulation snapshots reveal two com- 
mon structures, which are illustrated in Figure 01 In 
'pinned' structures (Figure H^), a single extended ICOF 
structure is folded on itself so that non-adjacent (non- 
intcrpcnctrating) icosahedra can be found to share ver- 
tices, edges, and faces with one another, pinning the 
fold, and producing a densely packed, stable, extended 
structure. A 'bridged' structure (Figure |3b) consists 
of proximate ICOFstructures that are joined by non- 
interpenetrating connections. The roles of these bridg- 
ing connections are of considerable interest since they 
serve to unify icosahedral networks and thus potentially 
affect the structure-property relationships of the MG. To 
demonstrate the significance of these bridging connec- 
tions, we measured the weighted connectivities of the 
icosahedral networks. 

Consider the icosahedra of a Cu-Zr MG to be the nodes 
of a graph, G, and their connections as the edges. That 
is, if two icosahedra share a vertex, edge, or face, or if 
they interpenetrate, then there is an edge between these 
two nodes. A connected subgraph, V, of G is defined to 
be a collection of nodes and their edges such that any 
two nodes in V are joined by a path of edges and nodes, 
and that no edge connects to a node outside of V. The 
size of V is defined to be the number of nodes in V. In 
this paper, an isolated node (one without edges) will be 
considered a connected subgraph of size 1. Our graph, 
G, is thus the union of all of its connected subgraphs. 
Let {Vi}™ =1 be the set of G's n connected subgraphs and 
i be the corresponding set of the sizes of the con- 
nected subgraphs. S max is the size of the largest sub- 
graph^) {Smax = max ({Si}i=i))- The weighted connec- 
tivity of G is then defined to be 

n „ 
k—1 max 

Thus a graph with Wq = 1 corresponds to a completely 
connected graph. A heuristic example of a graph, its con- 
nected subgraphs, and its weighted connectivity is pro- 
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FIG. 5. (Color Online) The average weighted connectivity 
for the ICOI-icosahedral network (blue circles) and the full 
icosahedral network (red squares), which includes all connec- 
tion types, for (a) Cu64Zr36 and (b) CusoZrso 



vided in Figure [BJ 

Figures and [3b show the weighted connectivities 
of both the ICOI-networks and the icosahedral networks 
that incorporate all connection types. It should imme- 
diately be pointed out that the decrease in weighted 
connectivity for high temperatures does not imply that 
the networks are highly connected prior to supercooling. 
Rather, the relatively low degree of icosahedral order- 
ing at these temperatures is responsible for the reduced 
measured connectivities. In both compositions of Cu-Zr, 
the non-interpenetrating connections that develop dur- 
ing supercooling dramatically enhance the unification of 
the networks. The effect of these connections is most 
obvious in Cu64Zr 36 MG (Fig. |5k); the weighted con- 
nectivity of the icosahedral network decreases rapidly as 
the alloy is cooled from 1100 K down to 800 K (T g ). In 
contrast to this, the weighted connectivity of the ICOI- 




FIG. 6. An example graph, G, containing two connected sub- 
graphs of sizes 4 and 3, respectively. The weighted connec- 
tivity of G is thus: Wa = \ + § 
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FIG. 7. (Color Online) (a) The average maximum icosahedron-to-icosahedron distance within a connected subgraph of the 
icosahedral network of Cu64Zr36. The simulation box employs periodic boundaries, thus seemingly disjoint atoms are attached 
to the connected subgraph across the box boundaries. As indicated, the horizontal dashed lines denote y/3l, where I is the 
respective simulation cube side length. The blue triangles correspond to the longest connected subgraph lengths when only 
considering ICOI for forming edges of the graph. The error bars indicate the standard deviation in the maximum subgraph 
length, [(b) - (e)] Illustrative simulation snapshots^ depicting the longest connected subgraph of the full icosahedral network 
at temperatures (b) 1200 K, (c) 1000 K, (d) 900 K, and (e) T g (800 K). The red spheres depict center Cu atoms of (0,0,12,0) 
clusters. The blue spheres depict both Cu and Zr atoms that serve strictly as vertices of the (0,0,12,0) clusters. The perspective 
in (e) was changed in order to show the extent of the percolating subgraph. 



network is not so significantly affected by the supercool- 
ing process. This suggests that non-interpenetrating con- 
nections are integral in linking the disjoint pieces of the 
ICOI-network as the alloy approaches its glass phase. In 
the case of Cu64Zr 36 , these bridging connections lead to 
an almost completely unified icosahedral network. The 
same network-unification was observed when the system 
size was increased to N — 10, 000 atoms. The weighted 
connectivity of Cu5oZr 50 's icosahedral network (Fig. [5b) 
exhibits a similar enhancement in comparison to its ICOI 
network, though network-unification is not observed in 
this alloy. 

These findings suggest that non-interpenetrating con- 
nections are an important facet of icosahedral networking 
in Cu-Zr MGs and perhaps in other MGs in which icosa- 
hedral ordering is prevalent. The heavy presence of these 
connections amidst the most highly-connected icosahedra 
suggests that they may assist in the dynamical slowing of 
these clusters. It is also reasonable to expect that their 
roles in enhancing the density and structural unification 
of the icosahedral network's MRO-structures contribute 
significantly to the structural-mechanical properties that 
were previously attributed only to ICOI-structure o 42 ' 46 . 
Pins and bridges may be particularly important for com- 
positions of Cu-Zr MG for which the development of 
ICOI is less extensive, and yet are good glass formers. 
As indicated earlier, only about 60% of the icosahedra 
in CusoZrso MG participate in interpenetrating connec- 
tions, but its icosahedral network is further developed via 



other connection types (Fig(3|). We are also interested 
in studying whether the icosahedral network's weighted 
connectivity is simply related to the population density 
of icosahedra, or if approaching a completely unified net- 
work reveals more complex features. 



VI. STATIC LENGTH SCALE 

Quenching a supercooled liquid into a glass produces 
major changes in the system's dynamics, yet there is 
no obvious long range structural change that accompa- 
nies this transition. Despite this, rigorous bounds have 
been proven to exist between length and time scales in 
a glassy system which mandate that a growing length 
scale must accompany the diverging relaxation times of 
glass transitions 3 - 3 .. Recent research efforts provide evi- 
dence for various growing correlation length s 18 ' 19 ' 21-23 ' 34 . 
However, there are findings that these scales (and in par- 
ticular static length scales- i.e., those that can be as- 
certained from static snapshots of the system) do not 
increase, when approaching the glass transition, rapidly 
enough, as suggested by various theories , in order to ac- 
count for the dramatic increase in the relaxation times^ 9 -. 
Static correlation lengths have been studied via 'point-to- 
set' correlation s 18 ' 24 ' 25 , three-point correlations 2 ^, pat- 
tern repetition size^l, and approaches employ graph the- 
oretical tools 2 ^ and shear rigidity penetration depths. 
Currently, most of the intuitive length scales, especially 
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FIG. 8. (Color Online) [(a) and (b)] The average percent composition of Cu present in: a single icosahedral cluster (including 
the center atom), the full icosahedral network, the liquid- like matrix (all atoms not involved in the icosahedral network) for 
(a) Cu64Zr36 and (b) CusnZrsn. [(c) and (d)] Display the average percent-abundance of Zr found in vertex, edge, face-sharing, 
and interpenetrating connections. The dotted lines indicate the percent-composition of the shell of an icosahedral cluster for 
(c) Cu64Zr36 and (d) CusoZrso. The error bars indicate the error in the mean values. 



static ones, do not increase as drastically as does the re- 
laxation time upon the glass transition. As glasses are 
rigid, a length characterizing the scale on which the su- 
percooled liquid can support shear— may be anticipated 
to naturally increase as the system "freezes" into a rigid 
amorphous solid. 

Motivated by these quests and, most notably, by 
the icosahedral network-unification that we observe in 
Cu64Zr36 during its glass transition, we introduce a 
static length scale, which interestingly exhibits dramatic 
growth upon quenching. Our length scale is the longest 
icosahedron-to-icosahedron distance within a connected 
subgraph of the full icosahedral network. It may be con- 
ceivable that the percolating icosahedral structure and 
the related increase in the interpenetrating icosahedral 
unit structure that we find might be related to an in- 
creased rigidity of the system as it rapidly cooled towards 
the low temperature glassy state. Connected subgraphs 
that cross the periodic boundaries are appropriately 'un- 
folded' across the boundaries before the distance is mea- 
sured, and the maximum length is averaged across con- 
secutive simulation snapshots. Figure [7J depicts the rapid 
growth in the length scale as Cu64Zr36 is cooled to T g . 
From 1600 K down to T gi the longest ICOI-chain length 
grows by a factor of 13. For the full icosahedral networks 
of the N = 2, 000 atom and N = 10, 000 atom systems, 
the longest connected subgraph lengths increase by fac- 
tors of 19 and 31, respectively. As expected from the 



percolating nature of the network, and as indicated by 
the horizontal dashed lines in Fig. [JJ the extent by which 
the length scale grows is apparently limited only by the 
simulation cell size. 



VII. COMPOSITIONAL INHOMOGENEITY 

Finally, we characterized the average chemical compo- 
sitions of a single icosahedral cluster (including its center 
atom), the icosahedral network, and the different con- 
nection types. The results are summarized in Figure [5] 
The icosahedral networks of both Cu-Zr alloys were found 
to be Zr-rich relative to their overall system composi- 
tions. This leaves the remaining liquidlike matrix of the 
systems to be Cu-rich, as shown in Figures [SJa, and [SJd. 
Figure displays the average percent abundance of Zr 
in each type of connection. Comparing these ratios to 
the average icosahedral-shcll composition (also displayed 
in Figure |5t) reveals that non-interpenetrating connec- 
tions do not appear to occur randomly between icosa- 
hedra, rather, they preferentially incorporate Zr-atoms. 
Because Zr has a larger atomic volume than Cu, these 
connections increase the density of the icosahedral net- 
work, while preserving the local composition of a single 
icosahedral cluster. 
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VIII. CONCLUSION 



In summary, we used a semi-empirical potential to con- 
duct MD simulations of CusoZrso and Cu64Zr36 across 
liquid and glass-phase temperatures. Our initial assess- 
ment of icosahedral ordering and networking within these 
systems provides additional support for the popular idea 
that interpenetrating connections serve to create robust, 
formative icosahedral structures. Non-interpenetrating 
connections are shown, however, to become increasingly 
prevalent as the alloys approach their glass phases. View- 
ing icosahedra as members of a network, interpenetrating 
icosahedra frequently participate in non-interpenetrating 
connections as well, and distinctly serve as the nodes 
of the network with the most connections. In particu- 
lar, non-interpenetrating 'pinning' and 'bridging' connec- 
tions are prominent amongst ICOI-structures, and may 
be of particular structural importance in Cu-Zr MGs in 
which ICOI are less extensive. To further characterize 
the structure of the icosahedral networks, we introduced 



the measure of weighted connectivity and used it to show 
that the network of Cu64Zr36 is almost completely uni- 
fied, in contrast to its ICOI-structures. Accordingly, the 
maximum connected-cluster length lends itself as a static 
length scale, whose rapid growth upon quenching is ap- 
parently limited only by the simulation cube size. Lastly, 
the vertex, edge, and face-sharing connections are shown 
to preferentially incorporate Zr atoms, amidst icosahe- 
dral networks that are Zr-rich compared to the system 
compositions. 
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